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Abstract 

The supersymmetric extension to the Standard Model offers a promising cold 
dark matter candidate, the lightest neutralino. I will review the phenomenological 
and cosmological constraints on the supersymmetric parameter space and discuss 
the prospects for the detection of this candidate in both accelerator and direct 
detection searches. 
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1 Introduction 


It is well known that supersymmetric models with conserved -R-parity contain one new 
stable particle which is a candidate for cold dark matter (CDM) [1], There are very strong 
constraints, however, forbidding the existence of stable or long lived particles which are 
not color and electrically neutral. Strong and electromagnetically interacting LSPs would 
become bound with normal matter forming anomalously heavy isotopes. Indeed, there are 
very strong upper limits on the abundances, relative to hydrogen, of nuclear isotopes [ 2 ], 
n/rtH ^ 1CT 15 to 10 ” 29 for 1 GeV < m < 1 TeV. A strongly interacting stable relic 
is expected to have an abundance n/riji < 1CT 10 with a higher abundance for charged 
particles. There are relatively few supersymmetric candidates which are not colored and 
are electrically neutral. The sneutrino [3, 4] is one possibility, but in the MSSM, it has 
been excluded as a dark matter candidate by direct [5] and indirect [ 6 ] searches. In fact, 
one can set an accelerator based limit on the sneutrino mass from neutrino counting, 
mo > 44.7 GeV [7]. In this case, the direct relic searches in underground low-background 
experiments require rrio > 20 TeV [5]. Another possibility is the gravitino which is 
probably the most difficult to exclude. I will concentrate on the remaining possibility in 
the MSSM, namely the neutralinos. 

There are four neutralinos, each of which is a linear combination of the R — — 1, 
neutral fermions [ 1 ]: the wino IV 3 , the partner of the 3rd component of the SU(2)l gauge 
boson; the bino, B , the partner of the U(l)y gauge boson; and the two neutral Higgsinos, 
H i and H 2 . In general, neutralinos can be expressed as a linear combination 


X = aB + /3W 3 + ytfi + SH-2 


( 1 ) 


The solution for the coefficients a,P, 7 and 5 for neutralinos that make up the LSP can 
be found by diagonalizing the mass matrix 
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where M l (M 2 ) are the soft supersymmetry breaking U(l) (SU(2)) gaugino mass terms. // 
is the supersymmetric Higgs mixing mass parameter and since there are two Higgs doublets 
in the MSSM, there are two vacuum expectation values, V\ and v 2 . One combination of 
these is related to the Z mass, and therefore is not a free parameter, while the other 
combination, the ratio of the two vevs, tan/3, is free. 

The most general version of the MSSM, despite its minimality in particles and interac¬ 
tions contains well over a hundred new parameters. The study of such a model would be 
untenable were it not for some (well motivated) assumptions. These have to do with the 
parameters associated with supersymmetry breaking. It is often assumed that, at some 
unification scale, all of the gaugino masses receive a common mass, m\/ 2 . The gaugino 
masses at the weak scale are determined by running a set of renormalization group equa¬ 
tions. Similarly, one often assumes that all scalars receive a common mass, mo, at the 
GUT scale (though one may wish to make an exception for the Higgs soft masses). These 
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too are run down to the weak scale. The remaining supersymmetry breaking parameters 
are the trilinear mass terms, A 0 , which I will also assume are unified at the GUT scale, 
and the bilinear mass term B. There are, in addition, two physical CP violating phases 
which will not be considered here. 

The natural boundary conditions at the GUT scale for the MSSM would include //, 
the two soft Higgs masses (mi and m 2 ) and B in addition to mp 2 , m 0 , and A 0 . In this 
case, by running the RGEs down to a low energy scale, one would predict the values 
of M z , tan/?, and the Higgs pseudoscalar mass, rrp (in addition to all of the sparticle 
masses). Since M z is known, it is more useful to analyze supersymmetric models where 
M z is input rather than output. It is also common to treat tan/? as an input parameter. 
This can be done at the expense of shifting /i (up to a sign) and B from inputs to outputs. 
When the supersymmetry breaking Higgs soft masses are also unified at the GUT scale 
(and take the common value mo), the model is often referred to as the constrained MSSM 
or CMSSM. Once these parameters are set, the entire spectrum of sparticle masses at the 
weak scale can be calculated. 


2 The Relic Density 

The relic abundance of LSP’s is determined by solving the Boltzmann equation for the 
LSP number density in an expanding Universe. The technique [8] used is similar to that 
for computing the relic abundance of massive neutrinos [9]. The relic density depends 
on additional parameters in the MSSM beyond mi/ 2 ,/q and tan/?. These include the 
sfermion masses, my, as well as tua, all derived from mo, A 0 , and mi/ 2 . To determine the 
relic- density it is necessary to obtain the general annihilation cross-section for neutralinos. 
In much of the parameter space of interest, the LSP is a bino and the annihilation pro¬ 
ceeds mainly through sfermion exchange. Because of the p-wave suppression associated 
with Majorana fermions, the s-wave part of the annihilation cross-section is suppressed by 
the outgoing fermion masses. This means that it is necessary to expand the cross-section 
to include p-wave corrections which can be expressed as a term proportional to the tem¬ 
perature if neutralinos are in equilibrium. LInless the neutralino mass happens to lie near 
near a pole, such as rn x ~ m z j 2 or m/,/2, in which case there are large contributions to 
the annihilation through direct s-channel resonance exchange, the dominant contribution 
to the BB annihilation cross section comes from crossed t-channel sfermion exchange. 

Annihilations in the early Universe continue until the annihilation rate T ~ &vn x 
drops below the expansion rate given by the Hubble parameter, H. For particles which 
annihilate through approximate weak scale interactions, this occurs when T ~ m x j 20. 
Subsequently, the relic- density of neutralinos is fixed relative to the number of relativistic 
particles. As noted above, the number density of neutralinos is tracked by a Boltzmann- 
like equation, 

- {(rv) (- n 2 - nl) (3) 

where n 0 is the equilibrium number density of neutralinos. By defining the quantity 
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( 4 ) 


/ = n/T 3 , we can rewrite this equation in terms of x — T/m x , as 
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where iV is the number of relativistic degrees of freedom. The solution to this equation at 
late times (small x) yields a constant value of /, so that n oc T 3 . The final relic density 
expressed as a fraction of the critical energy density can be written as [1] 


£l x h 2 ~ 1.9 
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where (T//T 7 ) 3 accounts for the subsequent reheating of the photon temperature with 
respect to y, due to the annihilations of particles with mass m < Xfm x [10]. The subscript 
/ refers to values at freeze-out, i.e., when annihilations cease. The coefficients a and b 
are related to the partial wave expansion of the cross-section, av — a + bx + .. .. Eq. (5 ) 
results in a very good approximation to the relic density expect near s-channel annihilation 
poles, thresholds and in regions where the LSP is nearly degenerate with the next lightest 
supersymmetric particle [11], 

The preferred range of the relic LSP density is provided by data on the cosmic mi¬ 
crowave background (CMB), which have recently been used to obtain the following 95% 
confidence range: Dcdm ^ 2 = 0.12 ± 0.04 [12], Values much smaller than Hcdm b 2 = 0.10 
seem to be disfavoured by earlier analyses of structure formation in the CDM framework, 
so we restrict our attention to n CDM h 2 >0.1. However, one should note that the LSP 
may not constitute all the CDM, in which case Dlsp could be reduced below this value. 
On the upper side, it is preferable to remain very conservative, in particular because the 
upper limit on Olsp sets the upper limit for the sparticle mass scale. Here, Ocdm ^ 2 < 0.3 
is used, while being aware that the lower part of this range currently appears the most 
plausible. 

Fig. 1 illustrates qualitatively the regions of the (mi/2,mo) plane favoured by LEP 
limits, particularly on b —» sy and cosmology. Electroweak symmetry breaking is 
not possible in the dark-shaded triangular region in the top left corner, and the lightest 
supersymmetric particle would be charged in the bottom right dark-shaded triangular 
region. The experimental constraints on m/> and b —» sy exert pressures from the left, 
as indicated, which depend on the value of tan /3 and the sign of ji. The indication of a 
deviation from the Standard Model in g^ — 2 disfavours /i < 0 and large values of m 0 and 
m 1/2 for fjL > 0. The region where Vt x h 2 falls within the preferred range is indicated in 
light shading, its exact shape being dependent on the value of tan f3. As discussed later 
in more detail, in addition to the ‘bulk’ region at low m o and mi/2, there are a number of 
regions at large values of mp2 and/or m 0 where the relic density is still compatible with 
the cosmological constraints. At large values of mm, the lighter stau, becomes nearly 
degenerate with the neutralino and coannihilations between these particles must be taken 
into account [13, 14], For non-zero values of Aq, there are new regions for which y — t 
coannihilations are important [15]. At large tan/3, as one increases mm, the pseudo¬ 
scalar mass, rriA begins to drop so that there is a wide funnel-like region (at all values of 
mo) such that 2 m x ~ rriA and s-channel annihilations become important [16, 17]. Finally, 
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Figure 1: Qualitative overview of a generic (mi/ 2 ,m 0 ) plane. The light (turquoise) shaded 
area is the cosmologically preferred region with 0.1 < Q x h 2 < 0.3. In the dark (brick red) 
shaded region at bottom right, the LSP is the charged f\, so this region is excluded. Elec- 
troweak symmetry breaking is not possible in the dark (pink) shaded region at top left. The 
LEP experimental constraints, in particular that on mh, and measurements ofb—ts'y ex¬ 
ert pressure from the left side. The BNL E821 measurement of g tl — 2 favours relatively 
low values of m 0 and mi/ 2 for p > 0. 


there is a region at very high mo where the value of p begins to fall and the LSP becomes 
more Higgsino-like. This is known as the ‘focus point’ region [18]. Also shown are the 
position of several benchmark points chosen for detailed phenomenological study [19]. 

Let us first focus on the ‘bulk’ region in the CMSSM for tan /3 = 10 and p > 0 
shown in Fig. 2 [13]. The light shaded region correspond to 0.1 < Q x h 2 < 0.3. The 
dark shaded region has m~ < m x and is excluded. The light dashed contours indicate 
the corresponding region in Pl x h 2 if one ignores the effect of coannihilations. Neglecting 
coannihilations, one would find an upper bound of ~ 450 GeV on mp 2 , corresponding to 
an upper bound of roughly 200 GeV on mp. 

Coannihilations are important when there are several particle species i, with different 
masses, and each with its own number density ni and equilibrium number density no,*. 
In this case [11], the rate equation (3) still applies, provided n is interpreted as the total 
number density, 

n = J2 n i , ( 6 ) 

i 

no as the total equilibrium number density, 

n 0 = J2 n o,i, (7) 

i 


4 













m 


1/2 


Figure 2: The light-shaded ‘bulk’ area is the cosmologically preferred region with 
0.1 < fl x h 2 < 0.3. The light dashed lines show the location of the cosmologically preferred 
region if one ignores coannihilations with the light sleptons. In the dark shaded region in 
the bottom right, the LSP is the f\, leading to an unacceptable abundance of charged dark 
matter. Also shown is the isomass contour m x ± = 104 GeV and m-h = 110,114 GeV, as 
well as an indication of the slepton bound from LEP. 


and the effective annihilation cross section as 


(o'eff^rel) 
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( 8 ) 


In eq. (4), m x is now understood as the mass of the lightest sparticle under consideration. 

Note that this implies that the ratio of relic densities computed with and without 
coannihilations is, roughly, 


R = 





( 9 ) 


where a = a+bx/2 and sub- and superscripts 0 denote quantities computed ignoring coan¬ 
nihilations. The ratio x^/x f ~ l+x° hi(g e ffCr e ff/(/i(Xo), where g e g = Yli 9 i{ m i/mi)‘ i ^ 2 e^ ( ' rni ~ ini ^ T . 
For the case of three degenerate slepton NLSPs, g e s = Hi9i — 8 and x®jx f ~ 1.2. 

The effect of coannihilations is to create an allowed band about 25-50 GeV wide in 
m 0 for my2 ^ 1400 GeV, which tracks above the rn^ = rn x contour. Along the line 
mf x = m x , R ~ 10, from (9) [13]. As mo increases, the mass difference increases and the 
slepton contribution to d e fr falls, and the relic density rises abruptly. 


3 The CMSSM with Universal Higgs masses 

A larger view of the tan /3 = 10 parameter plane is shown in Fig. 3 [21, 17, 19, 22, 23]. 
Included here are the most important phenomenological constraints (shown schematically 
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Figure 3: The (mi/ 2 ,mo) planes for (a) tan/? = 10 and /x < 0, (b) tan/? = 10 and ji > 0, 
assuming A 0 = 0 ,m t = 175 GeV and = 4.25 GeV. The near-vertical (red) 

dot-dashed lines are the contours rrifr = 114 GeV as calculated using FeynHiggs [20], 
and the near-vertical (black) dashed line is the contour m x ± = 103.5 GeV, shown only in 
(b). Also shown in (b) by the dot-dashed curve in the lower left is the corner excluded 
by the LEP bound of mg > 99 GeV. The medium (dark green) shaded region is excluded 
by b —> s'y, and the light (turquoise) shaded area is the cosmologically preferred regions 
with 0.1 < Vt x h 2 < 0.3. In the dark (brick red) shaded region, the LSP is the charged f\. 
The region allowed by the E821 measurement of a^ at the 2-a level, is shaded (pink) and 
bounded by solid black lines, with dashed lines indicating the 1-a ranges. 


in Figure 1). These include the constraints on the MSSM parameter space that are 
provided by direct sparticle searches at LEP, including that on the lightest chargino y ± : 
m x ± > 103.5 GeV [24], and that on the selectron e: mg > 99 GeV [25]. Another important 
constraint is provided by the LEP lower limit on the Higgs mass: rrin > 114.4 GeV 
[26] in the Standard Model. Since m./,, is sensitive to sparticle masses, particularly my, 
via loop corrections, the Higgs limit also imposes important constraints on the CMSSM 
parameters, principally m^. 

The constraint imposed by measurements of b —> sy [27] also excludes small values 
of mi/2- These measurements agree with the Standard Model, and therefore provide 
bounds on MSSM particles, such as the chargino and charged Higgs masses, in particular. 
Typically, the b —> sy constraint is more important for p < 0, but it is also relevant for 
H > 0, particularly when tan/? is large. 

The latest value of the anomalous magnetic moment of the muon reported [28] by the 
BNL E821 experiment is also taken into account. The world average of | (g /t — 2) 
now deviates by (33.9 ± 11.2) x 1CT 10 from the Standard Model calculation of [29] using 
e + e~ data, and by (17 ± 11) x 1CT 10 from the Standard Model calculation of [29] based 
on r decay data. Other recent analyses of the e + e~ data yield similar results. On some of 
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the subsequent plots, the formal 2-cr range 11.5 x 1CT 10 < 5a M < 56.3 x 1CT 10 is displayed. 
As one can see, the region preferred by g — 2 overlaps very nicely with the ‘bulk’ region 
for tan (3 = 10 and fi > 0. 


tan B = 35 , u. < 0 tan 3 = 50 , u. > 0 



m 1/2 (GeV) m 1/2 (GeV) 


Figure 4: As in Fig. 3 for (a) tan/1 = 35 and fi < 0 and (b) tan (3 = 50 and fi > 0 

As noted above, another mechanism for extending the allowed CMSSM region to 
large m x is rapid annihilation via a direct-channel pole when m x ~ \rriA [16, 17]. Since 
the heavy scalar and pseudoscalar Higgs masses decrease as tan f3 increases, eventually 
2 m x — rrtA yielding a ‘funnel’ extending to large and m 0 at large tan (3, as seen 
in Fig. 4. The difficulty and necessary care involved in calculations at large tan [3 were 
discussed in [17]. For related CMSSM calculations see [30]. 

Another way to see the effects of increasing tan [3 is plot look at the tan/3 — mi/ 2 
plane for fixed mo- This is shown in Fig. 5. As one can see, there is a charged A LSP 
in the upper right corner for both values of mo, though the area is reduced when mo is 
increased as is expected. One also sees the cosmological region moving up in tan f3 as m 0 
is increased. In b), one also sees the well defined coannihilation tail. Here, the importance 
of b —> s'y at large tan f3 is clearly displayed. 


4 The CMSSM with Non-Universal Higgs masses 

If one relaxes the unification condition for the two soft Higgs mass, m\ and m 2 , one could 
decide to make /i and rriA inputs and instead making the Higgs soft masses outputs. This 
is referred to as the CMSSM with non-universal Higgs masses or NUHM [31, 32, 23]. The 
attractive feature of radiative symmetry breaking in the CMSSM is maintained by using 
the electroweak minimization conditions to solve for mi and ?ri 2 , given values of /i and 
mA at the weak scale. 
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Figure 5: The tan (3 - mi/ 2 plane for (a) mo 
H > 0. All shading and line types are as in t 
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100 GeV and (b) m 0 = 200 GeV both for 
previous figures. 


The electroweak symmetry breaking conditions may be written in the form: 

m\{Q) = m\(Q) + m\{Q) + 2 p 2 (Q) + A A {Q) (10) 


and 


h 2 = 


m\ — ml tan 2 (3 + |m|(l — tan 2 (3) + Ajk 


( 11 ) 


tan 2 (3 — 1 + A^ 

where A^ and A^ 1,2 ) are loop corrections [33, 34, 35] and = mi^(m z ). The known 
radiative corrections [33, 36, 37] ci, c 2 and relating the values of the NUHM parameters 
at Q to their values at m z are incorporated: 


m i(Q) — rnj + Ci 

ml(Q) = m\ + c 2 

T 2 (Q) — [J^ + Cp. 


( 12 ) 


Solving for m\ and m 2 , one has 


m 2 (l + tan 2 (3) = m\{Q) tan 2 (3 — yU 2 (tan 2 (3 + 1 — A^) — (ci + c 2 + 2c M ) tan 2 f3 


-Aa(Q)taii 2 ,3 - im|(l - tail 2 (3) - A™ 


(13) 


and 


m 2 (l + tan 2 (3) = m 2 A (Q ) - /r 2 (tan 2 (3 + 1 + A| ( 2) ) - (d + c 2 + 2c M ) 

-A a (Q) + 1 “ tan 2 /?) + AW, 
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( 14 ) 
















which are used to perform the numerical calculations [23]. 

It can be seen from (13) and (14) that, if rriA is too small or /i is too large, then m 2 
and/or m| can become negative and large. This could lead to m\{Mx) + p 2 (Mx) < 0 
and/or m\(M x ) + /i 2 (M x ) < 0 , thus triggering electroweak symmetry breaking at the 
GUT scale. The requirement that electroweak symmetry breaking occurs far below the 
GUT scale forces one to impose the conditions m 2 (Mx) + p(M x ), m|(Mx) + p(M x ) > 0 
as extra constraints, which we call the GUT stability constraint 2 . 

The NUHM parameter space was recently analyzed [23] and a sample of the results 
found is shown in Fig. 6 . While much of the cosmologically preferred area with p < 0 is 
excluded, there is a significant enhancement in the allowed parameter space for p > 0 . 


tan (3 = 10, m 1/2 = 300 , m 0 = 100 



1000 


tan p = 10 , m 1/2 = 500 , m 0 = 100 



-2000 


2000 
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Figure 6: The NUHM (/qm^) planes for tan/? = 10, (a) mo = 100 GeV and m. 1/2 = 
300 GeV and (b) mo = 100 GeV and m. 1/2 = 500 GeV, with A 0 = 0. The (red) dot- 
dashed lines are the contours mh = 114 GeV, and the near-vertical (black) dashed lines 
are the contours m x ± = 103.5 GeV. The dark (black) dot-dashed lines indicate the GUT 
stability constraint. Only the areas inside these curves (small p) are allowed by this 
constraint. The light (turquoise) shaded areas are the cosmologically preferred regions 
with 0.1 < Vt x h 2 < 0.3. The dark (brick red) shaded regions is excluded because a charged 
particle is lighter than the neutralino, and the darker (dark blue) shaded regions is excluded 
because the LSP is a sneutrino. The medium (green) shaded region is excluded by b —> sy. 
The regions allowed by the E821 measurement of at the 2-a level are shaded (pink) and 
bounded by solid black lines. The solid (blue) curves correspond to m x = mA/2. 


As usual, there are dark (red) regions where there is one or more charged sparticle 
lighter than the neutralino y so that y is no longer the LSP. First, there are ‘shark’s 
teeth’ at \p\ ~ 300 GeV, m^ < 300 GeV in panel (a) of Fig. 6 where the f\ is the LSP. 
At small /i|, particularly at small m^ when the mass difference m 2 — m 2 is small, the tr 

2 For a different point of view, however, see [38]. 
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mass is driven small, making the f\ the LSP. At even smaller /i|, however, the lightest 
neutralino gets lighter again, since m x ~ // when // < M\ ~ Oimm. In addition, there 
are extended dark (red) shaded regions at large [//[ where left-handed sleptons become 
lighter than the neutralino. However, the electron sneutrino v e and the muon sneutrino z> M 
(which are degenerate within the approximations used here) have become joint LSPs at 
a slightly smaller |/t|. Since the possibility of sneutrino dark matter has been essentially 
excluded by a combination of V counting’ at LEP, which excludes mo < 44.7 GeV [7], 
and searches for cold dark matter, which exclude heavier z> weighing < 1 TeV [4], we 
still demand that the LSP be a neutralino x■ The darker (dark blue) shaded regions are 
where the sneutrinos are the LSPs, and therefore excluded. In the strips adjacent to the 
u efi LSP regions, neutralino-sneutrino coannihilation is important in suppressing the relic 
density to an acceptable level [23]. 

The thick cosmological region at smaller n corresponds to the ‘bulk’ region familiar 
from CMSSM studies. The two (black) crosses indicate the position of the CMSSM points 
for these input parameters. Extending upward in rriA from this region, there is another 
light (turquoise) shaded band at smaller /i|. Here, the neutralino gets more Higgsino- 
like and the annihilation to W + W~ becomes important, yielding a relic density in the 
allowed range 3 . For smaller /i|, the relic density becomes too small due to the X, ~ x' ~ 
X + coannihilations. For even smaller [^1 (< 30 GeV) many channels are kinematically 
unavailable and we are no longer near the h and Z pole. As a result the relic density may 
again come into the cosmologically preferred region. However, this region is excluded by 
the LEP limit on the chargino mass. 

The unshaded regions between the allowed bands have a relic density that is too high: 
Q x h 2 > 0.3. However, the f and v coannihilation and bulk bands are connected by 
horizontal bands of acceptable relic- density that are themselves separated by unshaded 
regions of low relic density, threaded by solid (blue) lines asymptoting to uia ~ 250 GeV. 
These lines correspond to cases when m x ~ m^/2, where direct-channel annihilation: 
X + y —> A, H is important, and suppresses the relic density [17, 16] creating ‘funnel’-like 
regions. 

In panel (b), a larger value of mi/ 2 is chosen. In this case, the previous region excluded 
by the neutral LSP constraint at large \fi\, migrates to larger |/t|. The ‘shark’s teeth’ for 
moderate \/i\ grow, reaching up to tua ~ 800 GeV. These arise when one combines a 
large value of with a relatively small value of mo, and one may find a f or even a 
e LSP. The large value of mi/ 2 also keeps the rate of b —> sy under control unless tua is 
small. The chargino constraint is similar to that in panel (a), whereas the mj t constraint 
is irrelevant due to the large value of m^. Finally, the GLIT mass-squared positivity 
constraint now allows larger values of \n\ < 1000 GeV. 

It is also interesting to examine the mo - m\/ 2 plane in the NLIHM [23]. In Fig. 7a, 
there is a bulk region at m\/ 2 ~ 50 GeV to 350 GeV, mo ~ 50 GeV to 150 GeV. As 
in the CMSSM, the f\ is the LSP in the bigger area at larger mm, and there are light 
(turquoise) shaded strips close to these forbidden regions where coannihilation suppresses 
the relic density sufficiently to be cosmologically acceptable. Further away from these 
regions, the relic density is generally too high. However, for larger m\/ 2 there is another 

3 This is similar to the focus-point region [18] in the CMSSM. 
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Figure 7: Compilations of phenomenological constraints on the MSSM with NUHM in 
the (mi/ 2 ,m 0 ) plane for tan/3 = 10 and (a) /i = 400 GeV, uia = 400 GeV, (h) /i = 
400 GeV, mA = 700 GeV, again assuming A 0 = 0. The shading is consistent with 
previous figures. The (blue) solid line is the contour m x = uia/^, near which rapid 
direct-channel annihilation suppresses the relic density. The dark (black) dot-dashed lines 
indicates when one or another Higgs mass-squared becomes negative at the GUT scale: 
only intermediate values of mi / 2 are allowed. 


suppression, discussed below, which makes the relic density too low. At small mi / 2 and 
mo the left handed slcptons, and also the sneutrinos, become lighter than the neutralino. 
The darker (dark blue) shaded area is where a sneutrino is the LSP. 

The near-vertical dark (black) dashed and light (red) dot-dashed lines in Fig. 7 are 
the LEP exclusion contours m x ± > 104 GeV and uih >114 GeV respectively. As in the 
CMSSM case, they exclude low values of mi/ 2 , and hence rule out rapid relic annihilation 
via direct-channel h and Z° poles. The solid lines curved around small values of mu 2 and 
m 0 bound the light (pink) shaded region favoured by and recent analyses of the e + e - 
data. 

A striking feature in Fig. 7(a) when mu 2 ~ 500 GeV is a strip with low kl x h 2 , which 
has bands with acceptable relic density on either side. The low -Q x h 2 strip is due to rapid 
annihilation via the direct-channel A, H poles which occur when m x = m^/2 = 200 GeV, 
indicated by the near-vertical solid (blue) line. The analogous rapid-annihilation strips 
occur in the CMSSM but at larger tan/3 as seen in Fig. 4. There, they are diagonal 
in the (m 1; / 2 ,mo) plane, reflecting a CMSSM link between mo and tua that is absent in 
this implementation of the NUHM. The right-hand band in Fig. 7(a) with acceptable 
Vt x li 2 is broadened because the neutralino acquires significant Higgsino content, and the 
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relic- density is suppressed by the increased W + W~ production. As mi/ 2 increases, the 
neutralino becomes almost degenerate with the second lightest neutralino and the lighter 
chargino, and the X~ x! ~ X ± coannihilation processes eventually push Q x h 2 < 0.1 when 
nr 1/2 > 700 GeV. 

The two dark (black) dash-dotted lines in Fig. 7(a) indicate where scalar squared 
masses become negative at the input GUT scale for one of the Higgs multiplets. One of 
these GUT stability lines is near-vertical at mm ~ 600 GeV, and the other is a curved 
line at ro 1 / 2 ~ 150 GeV, m 0 ~ 200 GeV. 

Panel (b) of Fig. 7 is for /t = 400 GeV and tua = 700 GeV. We notice immediately 
that the heavy Higgs pole and the right-hand boundary of the GUT stability region move 
out to larger mi/ 2 ~ 850,1050 GeV, respectively, as one would expect for larger tua- At 
this value of ttta, the right side of the rapid annihilation (‘funnel’) strip has disappeared, 
due to enhanced chargino-neutralino coannihilation effects. 

These two examples serve to demonstrate that the (roi/ 2 ,mo) plane may look rather 
different in the CMSSM from its appearance in the CMSSM for the same value of tan /?. 
In particular, the locations of rapid-annihilation funnels and X ~ X,' ~ X ± coannihilation 
regions are quite model-dependent, and the GUT stability requirement may exclude large 
parts of the (mi/ 2 , m 0 ) plane. 

4.1 Detection 

Because the LSP as dark matter is present locally, there are many avenues for pursuing 
dark matter detection. Direct detection techniques rely on an ample neutralino-nucleon 
scattering cross-section. The effective four-fermion lagrangian can be written as 

£ = Xl ll l b XQiln{.ocu + a 2 i7 5 )qi 
+ a 3iXXQiQi + a 4:iXl/ 5 XQilf 5 Qi 

+ a 5iXXQi'J 5 Qi + a&iXl b Xmi (15) 

However, the terms involving a 4i , a 4i , a 5i , and a 6i lead to velocity dependent elastic cross 
sections. The remaining terms are: the spin dependent coefficient, a 2i and the scalar 
coefficient a 3 j. Contributions to a? 2i are predominantly through light squark exchange. 
This is the dominant channel for binos. Scattering also occurs through Z exchange but this 
channel requires a strong Higgsino component. Contributions to aqj are also dominated 
by light squark exchange but Higgs exchange is non-negligible in most cases. 

Fig. 8 displays contours of the spin-independent cross section for the clastic scattering 
of the LSP x on protons in the mi/ 2 ,mo planes for (a) tan/? = 10, /t, < 0, (b) tan/? = 
10, n > 0 [39]. The double dot-dashed (orange) lines are contours of the spin-independent 
cross section, and the contours a S i = 10~ 9 pb in panel (a) and &si — 10~ 12 pb in panel 
(b) are indicated. The LEP lower limits on and m x ±, as well as the experimental 
measurement of b —> sy for fi < 0, tend to bound the cross sections from above, as 
discussed in more detail below. Generally speaking, the spin-independent cross section 
is relatively large in the ‘bulk’ region, but falls off in the coannihilation ‘tail’. Also, we 
note also that there is a strong cancellation in the spin-independent cross section when 
jj, < 0 [40, 41], as seen along strips in panel (a) of Fig. 8 where rn i/ 2 ~ 500 GeV. In the 
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Figure 8 : Spin-independent cross sections in the (mi/ 2 ,mo) planes for (a) tan (3 = 10 ,p< 
0, (b) tan (3 = 10, /i > 0. The double dot-dashed (orange) curves are contours of the spin- 
independent cross section, differing by factors of 10 (bolder) and interpolating factors of 
3 (finer - when shown). For example, in (b), the curves to the right of the one marked 
ICC 9 pb correspond to 3 x 10 -1 ° pb and 10 -1 ° pb. 

cancellation region, the cross section drops lower than 10 -14 pb. All these possibilities for 
suppressed spin-independent cross sections are disfavoured by the data on g fl — 2, which 
favour values of mi / 2 and mo that are not very large, as well as p > 0 , as seen in panel 
(b) of Fig. 8 . Thus g^ — 2 tends to provide a lower bound on the spin-independent cross 
section. 

Fig. 9(a) illustrates the effect on the cross sections of each of the principal phenomeno¬ 
logical constraints, for the particular case tan (3 = 10 p > 0. The solid (blue) lines mark 
the bounds on the cross sections allowed by the relic-density constraint 0.1 < ii x h 2 < 0.3 
alone. For any given value of mi/ 2 , only a restricted range of mo is allowed. Therefore, 
only a limited range of m o, and hence only a limited range for the cross section, is allowed 
for any given value of m x . The thicknesses of the allowed regions are due in part to 
the assumed uncertainties in the nuclear inputs. These have been discussed at length in 
[41, 40]. On the other hand, a broad range of rn x is allowed, when one takes into account 
the coannihilation ‘tail’ region at each tan/I and the rapid-annihilation ‘funnel’ regions 
for tan (3 = 35,50. The dashed (black) line displays the range allowed by the b — > sy 
constraint alone. In this case, a broader range of m 0 and hence the spin-independent 
cross section is possible for any given value of m x . The impact of the constraint due to 
m/i is shown by the dot-dashed (green) line. Comparing with the previous constraints, we 
see that a region at low m x is excluded by m^, strengthening significantly the previous 
upper limit on the spin-independent cross section. Finally, the dotted (red) lines in Fig. 9 
show the impact of the g^ — 2 constraint. This imposes an upper bound on mi / 2 and 
hence m x , and correspondingly a lower limit on the spin-independent cross section. 
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Figure 9: Allowed ranges of the cross sections for tan (3 = 10 (a) n > 0 for spin- 
independent elastic scattering. The solid (blue) lines indicate the relic density constraint, 
the dashed (black) lines the b —> s'y constraint, the dot-dashed (green) lines the irih con¬ 
straint, and the dotted (red) lines the g^ — 2 constraint. The shaded (pale blue) region is 
allowed by all the constraints, (b) The allowed ranges of the spin-independent cross sec¬ 
tion for j_i > 0. The darker solid (black) lines show the upper limits on the cross sections 
obtained from mh and b —> s'y, and (where applicable) the lighter solid (red) lines show 
the lower limits suggested by g^ — 2 and the dotted (green) lines the lower limits from the 
relic density. 


This analysis is extended in panel (b) of Fig. 9 to all the values 8 < tan/? < 55 and 
we find overall that [39] 


2 x 10- 10 pb < a S i < 6 x 10“ 8 pb, (16) 

2 x 1CT 7 pb < a SD <10" 5 pb, (17) 

for p > 0. (asD is the spin-dependent cross-section not shown in the figures presented 
here.) As we see in panel (b) of Fig. 9, mh provides the most important upper limit on the 
cross sections for tan /3 < 23, and b —> sy for larger tan (3, with g fJ — 2 always providing a 
more stringent lower limit than the relic-density constraint. The relic density constraint 
shown is evaluated at the endpoint of the coannihilation region. At large tan (3 , the Higgs 
funnels or the focus-point regions have not been considered, as their locations are very 
sensitive to input parameters and calculational details [42], 

The results from a CMSSM and MSSM analysis [40, 41] for tan/? = 3 and 10 are 
compared with the most recent CDMS [43] and Edelweiss [44] bounds in Fig. 10. These 
results have nearly entirely excluded the region purported by the DAMA [45] experiment. 
The CMSSM prediction [40] is shown by the dark shaded region, while the NUHM case 
[41] is shown by the larger lighter shaded region. 

I conclude by showing the prospects for direct detection for the benchmark points 
discussed above [46]. Fig. 11 shows rates for the elastic spin-independent scattering of 
supersymmetric relics, including the projected sensitivities for CDMS II [47] and CRESST 
[48] (solid) and GENIUS [49] (dashed). Also shown are the cross sections calculated in the 
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mx [GeV] 

Figure 10: Limits from the CDMS [43] and Edelweiss [44] experiments on the neutralino- 
proton elastic scattering cross section as a function of the neutralino mass. The Edelweiss 
limit is stronger at higher m x . These results nearly exclude the shaded region observed by 
DAM A [45]. The theoretical predictions lie at lower values of the cross section. 
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Figure 11: Elastic spin-independent scattering of supersymmetric relics on protons cal¬ 
culated in benchmark scenarios[46], compared with the projected sensitivities for CDMS 
II [47] and CRESST [ 48 ] (solid) mid GENIUS [49] (dashed). The predictions of our 
code (blue crosses) and Neutdriver [50] (red circles) for neutralino-nucleon scattering 
are compared. The labels correspond to the benchmark points as shown in Fig. 1. 
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proposed benchmark scenarios discussed in the previous section, which are considerably 
below the DAMA [45] range (10~ 5 — 10 -6 pb). Indirect searches for supersymmetric dark 
matter via the products of annihilations in the galactic halo or inside the Sun also have 
prospects in some of the benchmark scenarios [46]. 
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